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Abstract 



A minimal model is constructed for two-dimensional fracture propagation. The heterogeneous process 
zone is presumed to suppress stress relaxation rate, leading to non-quasistatic behavior. Using the Yoffe so- 
lution, I construct and solve a dynamical equation for the tip stress. I discuss a generic tip velocity response 
to local stress and find that noise-free propagation is either at steady state or oscillatory, depending only on 
one material parameter. Noise gives rise to intermittency and quasi-periodicity. The theory explains the ve- 
locity oscillations and the complicated behavior seen in polymeric and amorphous brittle materials. I suggest 
experimental verifications and new connections between velocity measurements and material properties. 
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The dynamics of cracks propagating in amorphous 
brittle media focused extensive study since the forties, 
mostly through quasi-static approaches and energetic ar- 
guments. In spite of recently renewed interest there are 
several fundamental issues that seem difficult to resolve 
in any simple way. For example: The limiting crack ve- 
locity, predicted to be the Rayleigh wave speed (RWS) 
in the bulk [1], is observed to be only about half of that; 
The mechanism for crack initiation and arrest is poorly 
understood; And the occurence of velocity oscillations 
[2] is still a puzzle. At the heart of the problem is the fact 
that the system's behavior depends on the lengthscale. 
While it is evident that the atomistic behavior differs 
from the continuous, it is this author's opinion that even 
on the continuum scale the physics near the tip is distinct 
from that far away and therefore should be treated dif- 
ferently. This may explain an apparent discrepancy: On 
the one hand, since the bulk shear wave speed (SWS) is 
higher than the crack velocity, it is clear that far from the 
crack quasi-static arguments should work well because 
the field relaxes to its static form, a a p = Kf a p{0)/^/2wr, 
sufficiently fast. Here a a p is the stress tensor, r is the 
distance from the crack tip, f a p depends only on the az- 
imuthal angle 9 and K is the (time-dependent) stress 
intensity factor. On the other hand, the inability of 
quasi-static theories to account for the above phenom- 
ena suggests that much of the behavior is determined by 
the local dynamics at the tip and hence that the propa- 
gation is a far-from-cquilibrium process, indescribable by 
approaches that appeal to energy balancing. In the two 
scale picture the noncquilibrium dynamics act in effect to 
dress the tip singularity as seen from afar. The matching 
of the near and far fields at the crossover scale then yields 
the far-away behavior of K. A reasonable guess would be 
that the crossover scale is of the order of the size of the 
processing zone (PZ) in front of the propagating crack. 
While the far quasi-static field is well understood within 
linear elasticity, there is little understanding of the short 
range physics although a few phenomenological dynamic 
equations have been advanced [3] [4] to explain the limit- 
ing tip velocity. 

It has been conjectured [5][6]that the reason for the 
complicated short-range behavior is the heterogeneous 
and fluid structure of the processing zone PZ. This con- 
jecture may be supported by observations of extremely 
slow relaxation rates of the stress at the tip after arrest 
[7] , rates that are an order of magnitude below expecta- 
tion had the relaxation taken place at the bulk speed of 
sound. This implies that the waves that re-establish the 
stress field in the PZ travel at a speed, c, that is much 
lower than the bulk SWS, probably due to scattering 
from microvoids. 

The model proposed here concerns the short-range 
dynamics and takes on board several ingredients: the low 
value of c in the PZ, the occurrence of different stresses 
for crack initiation and arrest [5] [8], and, based on exist- 
ing observations, an assumed velocity response to the tip 



stress. These suffice to construct and solve a dynamic 
equation. The explicit form of the velocity-stress rela- 
tion is not required for most of the results obtained here, 
only its qualitative behavior. The model leads to either 
a steady-state propagation at a limiting velocity or an 
oscillatory behavior, with the selection between the two 
modes depending on the location of the suppressed speed 
c on the velocity response function. Introduction of noise 
due to microvoid distribution is shown to give rise to an 
intermittent propagation that can turn into a quasiperi- 
odic behavior. 

Consider a line crack (not necessarily straight) in 
a two-dimensional material. The PZ in front of the tip 
is modelled as an effective continuous medium with a 
reduced SWS, c. The dependence of the crack dynam- 
ics on material properties enters through a velocity re- 
sponse function, v(cr) [9], where a is the local stress at 
the tip in the forward direction. As the crack prop- 
agates, the field near the tip adjusts to the changing 
boundary at a rate that corresponds to c. Observations 
that steady state propagation is at about half the bulk 
RWS, combined with the fact that c is much lower than 
the bulk (homogeneous) RWS, implies that the tip ve- 
locity can momentarily exceed the local value of c. This 
is a basic assumption in what follows. I comment that 
this does not violate the energy-balance which holds for 
scales away from the PZ because near the tip the dy- 
namic response, via), is swifter than the global energy 
equilibration process. The measured behavior of v(a) is 
hysteretic with two material-dependent thresholds: ah, 
above which propagation initiates, and 07 < ah, to which 
the stress has to drop for the crack to arrest [5] [8]. For 
a > ah the velocity is also known to increase very slowly 
with stress[10][ll]. Fig. 1 shows a qualitative form of 
v(a) that is consistent with experimental observations. 
This local non-monotonic response differs from that in 
[4]which depends on energy eqilibration far from the tip. 
Its locality allows to find the dynamics without further 
assumptions. A local two-branch velocity can be derived 
from atomistic models [12]. 

To derive the equation of motion of the tip, let us 
start from the Yoffc solution for the forward field of a 
propagating crack of length a [13], 



(C + a)/VC(C + 2a)" , 



(1) 



where £ is the distance from the tip and a^ is the tensile 
stress applied perpendicular to the propagation axis far 
away from the crack. In what follows the stress is mea- 
sured in units of a^ and a — > a/a^ is dimcnsionless and 
> 1. This solution assumes that the singularity of the 
field is always at the tip, which is consistent with a quasi- 
static picture. Consider, however, a situation wherein 
the dynamic response constrains the tip to overtake the 
density waves that adjust the field. In this situation the 
singularity in the stress field does not coincide with the 
location of the tip and the tip's stress drops to below 
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the static value. The difference between the static and 
dynamic stresses at the tip depends on the tip's veloc- 
ity, v = dl/dt and the propagation history. The dynamic 
stress is found from (1) by putting £ = (I — ct)Q(l — ct), 
where I is the tip's position and the step-function, O, en- 
sures that when the shear wave catches up the tip stress 
stays at the static value. When O = the tip stress di- 
verges as expected and traditional quasi-static solutions 
apply [14]. Focusing on non-quasi-static propagation, I 
assume = 1 during the entire growth. When © alter- 
nates between and 1 one simply pieces the solutions 
together. Taking the time derivative of (1) we have 
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o = -Ca 2 /[C(C + 2a)] 
Using (1), we can invert relation (2): 

C/a=-a/(^-lf 2 



(2) 



(3) 



We define Q = (v — c) as ( = cu, where u is a reduced 
velocity. Upon substitution in (3) we can readily solve 
the equation: 



t - t = — 



a(t) 



ds 
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The kinetics are thus determined by the response func- 
tion through the stress-dependence of u(a). Relation (4) 
is the bare result of this report. It is an exact derivation 
from the Yoffe solution. It gives the general time depen- 
dence of the stress at the crack tip. Once the stress his- 
tory is found from this relation one substitutes it in v(&) 
to obtain the velocity history. We now proceed to analyse 
the consequences of this result, assuming the qualitative 
response shown in Fig. 1. It is convenient to classify 
the behavior in terms of the ratio A = cjv\. The rea- 
son is that, as is shown below, the mode of propagation 
depends only on this ratio. 

A > 1: The point (a(c), c) is on the upper branch 
of v{a). Suppose that initially a < 07. The velocity is 
momentarily zero (or very low) and the tip stress builds 
up to 07,. At this stage the system 'jumps' to the upper 
branch and fast motion ensues. From relation (4) we 
see that for a^ > a(c) (< <r(c)) the stress will decrease 
(increase) until a converges to a(c) whereafter the tip 
propagates at a velocity c and a fixed distance ahead of 
the density waves. Thus (<x(c), c) is a fixed point of the 
equation of motion. The behavior at the vicinity of this 
point can be found by linearization of relation (4): 
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(5) 

Since v(a) near a(c) is smooth and positive 7 is regu- 
lar and positive and the fixed point is stable, namely, 
a steady state propagation at a limiting velocity c is a 



stable fixed point of the dynamics. A typical such his- 
tory of v is shown in fig. 2. An interesting implication 
of this result is that the experimentally observed limit- 
ing crack velocities give in fact the value of c and hence 
the local stress relaxation rate. This suggests a check 
of this model by comparing the limiting velocity to the 
speed of sound in the PZ. It is intriguing to note that 
even in the absence of a global energy balance criterion 
the crack velocity converges to the SWS, albeit the local 
value, c. Observations that v increases very slowly with 
K [15] in this regime indicate a small value of dv/da 
along the upper branch. In view of the present analysis, 
this agrees with the reported velocity behavior immedi- 
ately after crack initiation [5] [8] [10]. Another check of 
this picture can be suggested: In some experiments a 
drop in the stress has been measured after crack initia- 
tion [16]. This suggests that in those systems a(c) < cr^, 
a conclusion that can be checked by independent meth- 
ods as a test of this analysis. 

A < 1: To analyse this case, let us assume again 
that initially the tip stress is lower than 07 . From relation 
(4) the stress will increase until it reaches ah, whereupon 
the crack will start propagating as for A > 1. The velocity 
and the stress will then gradually decrease. Since c < vi 
is not a point on the upper branch the system cannot 
settle into a steady state as before and at 07 it flips back 
to the lower branch. There the crack halts momentarily, 
the stress at the tip builds up again to ah and the cycle 
repeats itself. This is a relaxation cycle whose period is 
found from (4): 



if 



l/u ub (s) - l/ui b (s) 

(S 2 - l) 3 / 2 



ds 



(6) 



where u u b > and uu, < are, respectively, the values of 
u along the upper and lower branches. When the velocity 
vanishes along the lower branch uy, = — 1. A typical 
velocity history in this case is also shown in fig. 2. 

A = 1: This marginal case is sensitive to the value 
of dv/da at 07. If the derivative is regular one can easily 
see that the analysis is the same as for A > 1. The only 
difference is that a can only approach cr(c) = 07 from 
above because for a < a\ the only motion is up the lower 
branch. If dv/da diverges at 07 the behavior depends 
on the detailed form of the divergence. For illustration, 
consider the form 



v = vi exp 



, a 
a I 1 

07 



u ~ Const. + (a — aif 



(7) 

with < v < 1. The behavior near the fixed point is 
found by using (7) in Eq. (4), 



a(c) ~ (r - t) 1/(1 -^ 



(ro - t) 



u/(l- V ) 



(8) 

where t > t is a constant. Now the propagation rate 
converges to c as a power law rather than exponentially. 
It should be noted that this propagation mode is sensitive 
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to small fluctuations that can easily flip the system to the 
lower branch. A small noise in this case will give rise to 
a quasi-periodic motion similar to that discussed below. 

Intermittency and quasi-periodicity: Since the 
PZ is inhomogeneous one expects fluctuations in the lo- 
cal properties which may well go beyond the effective 
medium assumption. For simplicity, let us restrict the 
discussion only to fluctuations in the tip stress during 
propagation, a = ao(t)—T](t), where er is the stress in the 
absence of noise. This corresponds to a situation where 
the crack encounters microvoids of varying sizes along its 
path. Upon association of a microvoid to the crack the 
tip stress drops momentarily with the drop depending on 
the microvoid's size. The fluctuations in microvoid sizes 
give rise then to noise in the tip stress. A fluctuation 
during steady-state propagation that reduces a by more 
than S = <j{c) — ai (see fig. 1) flips the system to the lower 
branch. The system then has to go through the process 
of stress increase to <7h, jump to the upper branch and 
convergence to a(c) again. The time that this process 
takes depends on the original fluctuation, f] > 5, and can 
be found by applying (4) to the motion along the two 
branches: 
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(9) 

It is the occurrence frequency of the flips between the 
branches which determines to a large extent the observ- 
able behavior. This frequency depends both on the noise 
characteristics and the value of 5. The stochastic velocity 
behavior can be obtained from the statistics of r\ by using 
relation (9). For example, the probability density of T, 
P(T), can be found from that of n, Po(rj), by inverting 
relation (9) to obtain rj(T) and substituting in 



P(T) = P (v(T)) (d V /dT) . 



(10) 



From (9) one can also find the effects of various forms 
of the noise temporal correlations, < rj(t)q(t') > on the 
velocity history. A detailed analysis of the statistics, in- 
cluding the explicit dependence on the distribution of 
microvoid sizes, will be reported shortly. Here I only 
point out a few intriguing consequences for A > 1 and 
a{c)<Uh'- First, a low occurrence frequency (i.e., <l/r) 
of rj > 5 leads to an intermittent behavior wherein the tip 
is 'knocked' occasionally from the steady state and then 
returns to it only to be knocked out of it again at a later 
time. A plot of such a history is shown in fig. 3. Second, 
a very high occurence frequency of fluctuations rj > a^—ai 
gives rise to a quasi-periodic behavior as follows: As the 
tip stress builds up along the lower branch to a h the 
system flips to the upper branch. A fluctuation then im- 
mediately knocks the system back to the lower branch, 
not allowing it to settle into the steady state. The mean 
period will be then close to the time spent on the lower 
branch, namely, w (a/c) [(1- 1/af)- 1 / 2 - (1 - l/a^)" 1 / 2 ] . 



Thus the microvoid size distribution determines the ob- 
servable behavior by governing the statistics of 77. Since 
this distribution also plays a major role in determining 
the roughness of the fracture surfaces my analysis can 
rigorously link roughness measurements to the velocity 
history. 

To summarize, a minimal theoretical model has been 
proposed to explain the rich behavior observed in crack 
propagation in amorphous and polymeric materials. The 
theory is a direct consequence of the observations of slow 
relaxation rates of the tip stress, the occurence of dif- 
ferent initiation and arrest stresses, and the deduced 
qualitative form of v(a). The mode of propagation has 
been found to depend only on one material parameter, 
A = c/vi. For A > 1 the propagation speed saturates 
to a limiting value, while for A < 1 it oscillates pe- 
riodically. Noise gives rise to a spectrum of behavior 
ranging from intermittent to quasi-periodic propagation. 
The model explains naturally recent observations of os- 
cillations in polymeric materials and measurements for 
its validation have been suggested. It predicts that the 
low-noise steady-state growth rate is exactly c, the wave 
speed in the PZ, and should be possible to test experi- 
mentally. Low A is expected to correspond to high disor- 
der and vice versa. So by manipulating the disorder one 
may tune A. I should remark that the effective contin- 
uum approximation of the PZ probably breaks down for 
too broad a distribution of microvoid sizes, and a statis- 
tical treatment is more adequate. E.g., for propagation 
velocities of order c <~ 500m/s and given the fact that 
currently observations are limited to times of fisec and 
higher, the effective continuum assumption should hold 
when microvoids are smaller than 500/jm. Microvoids do 
not usually reach such sizes in polymeric materials and 
therefore this model should do a good job explaining ex- 
perimental observations [2] in these systems. A comple- 
mentary statistical analysis for broad microvoid distribu- 
tions is currently under way and will be reported shortly. 
Finally, many ramifications of this model remain to be 
explored: the effects of noise correlations < rj(t)n(t') > 
on the dynamics, the effects of realistic distribution of 
void sizes, and the implications of the statistics and ve- 
locity history on the surface roughness, to name a few. 

FIGURE CAPTIONS 

1. A generic plot of v(a). 

2. Typical velocity histories in the steady-state (solid) 
and periodic (dashed) regimes. 

3. A typical velocity history in the intermittent regime. 
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